AIM  4C2 


unclassified 


INVESTIGATION  OF  LASEA-INOUCED  PLASNA  PROCESSESIU) 
VISIDVNE  INC  BURLINGTON  MA  P  C  IP  ET  AL  IB  DEC  8£ 
AFGL-TR-87-BaS2  FlSS28-BT-C-ai49 


AD- A 190  463 


one  FILE  COPS 


AFGL-TR-87-0052 

INVESTIGATION  OF  LASER- INDUCED  PLASMA  PROCESSES 


Precila  C.F.  Ip 
Russell  A.  Armstrong 

Mission  Research  Corporation 
One  Tara  Blvd,  Suite  302 
Nashua  N.H.  03062 

and 

James  C.  Baird 

Brown  University 
Department  of  Physics 
Providence,  R.I.  02912 


Prepared  for; 

Visidyne  Corporation^^ 
10  Corporate  Place 
South  Bedford  St. 
Burlington,  MA  01803 


30  December  1986 


Scientific  Report  No.  14 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED 


AIR  EORCE  GEOPHYSICS  LABORATORY 
AIR  FORCE  SYSTEMS  COMMAND 
UNITED  STATES  AIR  EORCE 
HANSCOM  AFB,  MA  01731 


o 


This  technical  report  has  been  reviewed  and  is  approved  ior  publication 


UKAN  I'.  KIMUAl.l, 

Contract  i*laiia{;cr' 

Atiiioaphci'i  c  Backy rouiidb  Br’anch 
Itifi'at'cd  Tccliiuilugy  Division 


UOBEUT  k.  O'NEIL,  Chid' 
Infrared  Dynamics  lii'aiicli 
Infrared  Tcchnoloi'v  Division 


FOR  THE  COMMANDER 


KAVlliALL  E“  'HURPHY,  DirectoK  “ 
Infrared  Technology  Division 


This  report  has  been  reviewed  by  the  ESD  Public  Affairs  Office  (PA)  and  is 
releasable  to  the  National  Technical  Information  Service  (NTIS). 


Qualified  requestors  may  obtain  additional  copies  from  the  Defense 
Technical  Information  Center.  All  others  should  apply  to  die  National 
Technical  Information  Service. 


If  your  address  has  changed,  or  if  you  wish  to  be  removed  from  the  mailing 
list,  or  if  the  addressee  is  no  longer  employed  by  your  organization, 
please  notify  AFGL/DAA,  Hanscom  AFB,  MA  01731.  This  will  assist  us  in 
maintaining  a  current  mailing  list. 


Do  not  return  copies  of  this  report  unless  contractual  obligations  or 
notices  on  a  specific  document  requires  that  it  be  returned. 


UNCLASSIFIED 


REPORT  DOCUMENTATION  PAGE 


lC^S5>?'Ci’*'ON  DO'.VNGPAOiNG  SC:^6DL/L5 


3  ^lSTalBlj^0^4/ AVAiLAaiLi"'^  OF  REPORT 

Approved  for  public  release,  distribution 
uni inited 


■■  =”i:ORMiMG  ORGANIZAT  ON  REPORT  \lo.'V13fcR(S) 

MRC-NSH-R-0002 


5  .VtONITORING  ORGANIZA'CM  REPORT  'jlV13£R(S) 

AFGL-TR-87-0052 


Ea  '.A, ME  Dr  PERFORMING  ORGANIZATION  6b  OFFICE  SYMBOL  7a  NAME  OF  MONITORING  ORGANIZATION 
'.M  '  ■.  Rf  SEARCH  r.ORPORAT  '  ON  (if  ippheiOl*)  cnorc  nmouvc  t  rc  i  Anno  A' 


AIR  FORCE  GEOPHYSICS  LABORATORY 


“<-  iDGRSSS  Cry,  State  ara  ^IPCjdel 

;  TARA  BE  VI),  SI'ITE  302 
NASHUA,  NH  03002 


7b  ADDRESS -C/tv,  Sfare  and  ^IP  Code) 

HANSCOM  AFB,  MA  01731 


ia.  NAME  OF  FUNDING  -  SPONSORING 
DRGANiZAT'ON 


VISIDYNE  CORPORATICN 


3b  OFFICE  SYMBOL  I  9  PROCUREMENT  NSTRUMENT  OENTIFiCAT'ON  NUMBER 
(If  applicable)  I 

F19628-83-C-0149 


ic  ADDRESS  iCry.  State,  and  ZIPCode) 

10  CORPORATE  PLACE 
tOUTH  BEDFORD  STREET 
nilUl  IMI'TllM  MA  ni 


'  ~.z  'Jnctuae  Security  Claisification) 


1  to  SOURCE  OF  funding  numbers 

PROGRAM 

pROjEC' 

'ASK 

'.VORK  uNiT 

ELEMENT  NO 

NO 

NO 

ACCESSION  NO 

62101E 

7670 

10 

A.L 

INVESTIGATION  OF  LASER- INDUCED  PLASMA  PROCESSES 


'Z  personal  Aw'aOR(S) 

Dr.  Precila  C.F.  Ip,  Dr.  Russell  A.  Armstrong  and  Dr.  Dames  C.  Baird* 


■3a  '■■'Sf  CF  REPORT 

Scientific  No.  14 


'6  supplementary  notat.on 

*  Brown  University 


13b  time  COVERED 

14  DATE  OF  REPORT  .Year  Montt^,  Day) 

FROM  3-14-86  To1 

1986  December  30 

COSATi  CCOES 


'8  SUBJECT  ■^ERMS  (Continue  an  reverse  if  necessary  and  identify  by  block  number) 

GROUP  I  SUB-GROUP  I  Emissions,  Oxygen  Plasma,  Rydberg  States,  Nuclear  Effects 


'9  A8STRAC'  Continue  an  reverse  if  necessary  and  identify  by  block  number) 


Prelirinary  results  on  the  hydrodynamics  and  Stark  calculations  of  the  AFGL-LINUS  plasma  are 
prison :.?d.  Also  reported  is  a  kinetic  analysis  based  on  the  most  recent  spectroscopic 
tuiipcTctures .  Various  calculations  of  different  aspects  of  the  LINUS  plasma  are  shown. 
Fmaliv.  the  codes  (EXCEDE/ARTIC)  which  have  potential  applications  to  LABCEDE  are  discussed. 


;0  DlSTO'BU'  ON  .  availability  OF  A8STRAC' 

C -NCLAS3.F’ED/UNL  MITFD  21  SAME  as  RPt  □  Dy  C  USERS 

21  A3STRAC'  SECURITV  C^SSiF'CATiQN 

UNCLASSIFIED 

;:a  NAME  OF  RESPONSIBLE  NOiV'OUAL 

Dean  F.  Kimball 

22b  TeuEP'-'ONE  (Include  Area  Code)  | 

(617)  377-3642 

22c  OFF'CE  S-yMSOl 

AFGL/LSP 

DO  FORM  1473,  34  MAR 


33  APR  edition  may  be  used  until  eKnausted 
All  other  editions  are  obioiet#. 


SECURITY  CLASSiFiCA'iQN  QF  'UlS  PAGE 

UNCLASSIFIED 


L>i 


TABLE  OF  CONTENTS 


tg 

k  '^m 

•. 

R:-- 


1 1 1 


SECTION 

1 

INTRODUCTION 

PAGE 

1 

■ 

II 

HYDRODYNAMICS 

2 

1 

III 

STARK  CALCULATION 

31 

0  ,  ’• 

IV 

KINETIC  EFFECTS 

51 

■.■■- 

V 

MISCELLANEOUS  ISSUES 

58 

m 

VI 

EXCEDE  AND  ARCTIC  CODES 

72 

r."j 

'*k  ■ 

r/-' 

VII 

REFERENCES 

77 

I 


LIST  OF  FIGURES 

PAGE 


Mg.  1.  SALE  Calculation  for  the  Sample  Case  Running  on  the  CRAY.  4 

A.  Pressure  B.  Density  C.  Energy  D.  Speed. 

Mg.  2.  SALE  Calculation  for  the  Sample  Case  Running  on  the  PC.  8 

A.  Pressure  B.  Density  C.  Energy  D.  Speed. 

Mg.  3.  The  Computing  Mesh  being  Input  into  the  SALE  Code.  13 

Mg.  4.  SALE  Calculation  for  the  LINUS  Case  at  1  nsec.  14 

A.  Pressure  B.  Density  C.  Energy  D.  Speed. 

Fig.  5.  SALE  Calculation  for  the  LINUS  Case  at  3  nsec.  18 

A.  Pressure  B.  Density  C.  Energy  D.  Speed. 

Fig.  6.  SALE  Calculation  for  the  LINUS  Case  at  5  nsec.  22 

A.  Pressure  B.  Density  C.  Energy  D.  Speed. 

Fig.  7.  Bow  Tie  Geometry  for  Laser  Focal  Region  Illustrating  30 

Shock  Wave  Interaction. 

Fig.  8.  Zero  Plasma  Field  Spectra  for  18  Basis  States  33 

(4d;  4f;  5d;  5f;  5g;  6d;  6f;  6g;  6h;  7f;  7g;  7h;  7i;  8f; 

8g;  8h;  8i ;  8j) . 

Fig.  9.  01  Spectra  for  18  Basis  States  for  Electron  Densities  of  39 

A.  10^^,  B.  10^^,  and  C.  10^^  el ectrons/cm^ . 

Fig.  10.  Assumed  Temperature  Profile  for  Kinetic  Calculation.  54 

Fig.  11.  Calculated  log  density  of  species  as  a  function  of  time  55 

for  10  torr  oxygen. 

Fig.  12.  Calculated  log  density  of  species  as  a  function  of  time  56 

for  100  torr  oxygen. 

Fig.  13.  Calculation  of  Emission  Coefficients  for  Free-Free  63 

Transition  in  the  Field  of  an  Ion. 

Fig.  14.  Front  View  of  CVF.  65 

Fig.  15.  Wavelength  Shift  of  CVF  as  a  Function  of  Incident  Angle.  66 

Fig.  16.  Wavelength  Shift  of  CVF  as  a  Function  of  Temperature.  66 

Fig.  17.  Optical  System  Using  One  Lens  for  Using  CVF  in  LINUS.  68 

Fig.  18.  Optical  System  Using  Lens,  Collimator,  and  Slit  for  68 

Using  CVF  in  LINUS. 


i  V 


FIGURESfContinued^ 


Fig.  19.  EXCEDE  Calculation  for  the  Case  of  4.5  eV  Incident 
Electrons  and  10  Collisions  in  N2. 

Fig.  20.  Radial  Electron  Density  Profiles  in  LABCEDE  for  4.5 
Incident  Electrons  in  N2. 


.Acce^i'jlon  For 

NT  1 3  GK/u'cI 

1  DTIC  TAB 
j  Unarmciinced 
j  JuGtiflcatl 

□ 

!  By 

j  Distribution/ 

!  Aval 

List 

tl\ 

lahllll 

Avail 

Spec 

-y  Codas 

and/or 

lal 

V 


eV 


TABLES 


PAGE 

Table  1.  Comparison  of  Shock  Front  Velocity  Calculated  from  SALE  26 
and  from  Blast  Wave  Theory. 

Table  2.  Reaction  Set  for  LINUS  Recombination 
Table  3.  Spectroscopic  Temperatures  for  Oxygen  in  LINUS 


52 

53 


1.  INTRODUCTION 


v. 

L*. 


« 


This  final  report  summarizes  work  performed  during  the  period  March 
1986-December  1986  under  Visidyne  contract  No.  F19628-83-C-0149,  MRC  work 
order  15105  WPR. 

The  purpose  of  the  study  is  to  investigate  properties  of  the  laser- 
induced  plasma  in  the  LINUS  experiment  in  order  to  test  its  validity  as  a 
benchmark  source  for  a  plasma  plume  model  and  as  a  basis  for  more 
sophisticated  model  development.  Plasma  recombination  processes  are 
responsible  for  the  large  amount  of  radiation  (ultraviolet,  visible,  and 
infrared)  generated  in  a  plasma  initiated  by  a  laser  beam.  The  spectral 
information  gained  and  the  processes  responsible  for  their  generation  are 
important  in  their  potential  application  to  understanding  the  spectral 
emission  arising  from  an  actual  nuclear  airburst.  The  results  may  be 
benchmarked  against  Sappenfield's  plasma  radiance  model [1]  which  has  been 
extrapolated  to  the  infrared  regime  where  there  is  little  data  available. 
However,  before  any  comparison  can  be  made,  we  must  first  complete  a  full 
detailed  analysis  of  the  laser-induced  plasma. 

Section  fl  discusses  the  hydrodynamics  in  LINUS.  Section  III 
presents  a  Stark  calculation  for  application  to  analysis  of  LINUS. 
Section  IV  deals  with  the  kinetic  effects  and  Section  V  shows  the 
calculations  of  various  important  issues  in  LINUS. 

Section  VI  discusses  the  Excede  and  Arctic  codes  which  are  important 
in  their  application  to  LABCEDE,  a  LINUS-rel ated  low  density  plasma 
project.  The  LABCEDE  experiment  employs  an  electron  beam  to  excite 
atmospheric  gases  the  purpose  of  which  is  to  validate  and  develop 
atmospheric  weapons  effects  codes.  The  Arctic  code  calculates  radiance 
arising  from  electron  excitation  of  atmospheric  gases  in  an  aurora  and  may 
be  useful  in  the  analysis  of  linear  and  non-linear  effects  in  LABCEDE. 
The  primary  interest  in  relating  LINUS  and  LABCEDE  is  that  atomic  oxygen 
emissions  seen  in  LABCEDE  have  some  common  characteristics  with  those  in 
LINUS[2] . 
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II.  HYDRODYNAMICS  IN  LINUS 


A.  INTRODUCTION 

This  section  summarizes  the  preliminary  investigation  into  the 
dynamics  of  the  discharge  region  in  the  LINUS  experiment.  The 
hydrodynamics  study  is  important  for  understanding  the  general  morphology 
of  the  plasma  expansion  and  the  effects  on  recombination  and  radiative 
properties.  Furthermore,  gradient  profiles  of  electron  density  and 
tr-',  mature  are  required  to  characterize  the  laser  induced  plasma.  The 
important  hydrodynamical  processes  to  be  examined  are  the  generation  of 
shock  fronts  and  initial  expansion  due  to  laser  energy  deposition,  shock 
front  effects,  post-shock  cooling,  and  pressure  recovery. 

We  have  initiated  the  calculations  of  the  shock  front  and  expansion 
characteristics  assuming  simple  ordinary  fluid  dynamics.  The  code  which 
is  used  for  this  purpose  is  the  Simplified  Arbitrary  Lagrangian  Eulerian 
(SALE)  code  originally  written  by  Amsden  et  al.[3]  from  Los  Alamos 
Scientific  Lab  (now  LANL).  SALE  is  a  general  fluid  code  which  uses  a 
simplified  numerical  fluid  dynamics  computing  technique  to  calculate  two- 
dimensional  fluid  flows.  It  can  handle  flow  speed  from  the  incompressible 
to  the  supersonic  limit  in  the  Lagrangian  or  the  Eulerian  mode,  or 
anywhere  in  between.  In  the  Lagrangian  mode,  the  cell  vertices  move  with 
the  fluid  whereas  in  the  Eulerian  mode,  the  grid  vertices  remain 
stationary  with  respect  to  the  fluid.  The  code  has  been  widely  used  and 
well  tested.  Its  major  limitations  for  the  present  application  are  the 
treatment  of  only  one  species  and  no  allowance  (yet)  for  energy  loss. 

The  code  is  based  on  the  solution  to  the  Navier-Stokes  fluid  equation 
of  motion  and  the  mass  and  internal  energy  equations.  One  begins  by 
specifying  the  geometry  of  the  computing  mesh  which  is  a  set  of 
quadr i 1 ateral  cells.  For  each  cell,  the  code  calculates  momentum, 
position,  volume,  energy,  mass  density,  and  pressure.  The  pressure  is 
calculated  from  the  specified  equation  of  state.  Input  variables  are  the 
initial  time  step,  the  initial  density  and  specific  internal  energy.  All 
subsequent  time  steps  are  automatically  adjusted.  The  code  calculates  all 
the  above  parameters  for  every  cell  and  each  time  step  until  it  comes  to 
the  terminating  time  (input  by  user). 
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B.  RUNNING  ON  THE  PC 


The  code  resides  on  the  Los  Alamos  Cray  computer.  We  have  been  able 
to  convert  the  code  to  run  on  the  PC.  The  Cray  is  about  160  times  faster 
than  the  PC  but  to  run  a  job  is  considerably  more  expensive  and  less 
accessible.  The  SALE  code  is  about  2000  lines  and  the  executable  image  is 
about  340K  bytes  which  is  well  within  the  memory  buffer  {640K)  of  the  PC. 

In  the  test  case  the  initial  parameters  are  -1  /xJ  of  energy  deposited 
into  a  focal  region  of  50  /xm  by  120  /xm.  Although  the  initial  test  energy 
was  lower  than  that  in  LINUS,  the  morphological  results  should  be  similar. 
This  sample  case  was  run  on  both  the  CRAY  and  the  PC.  The  results  are 
identical  except  for  differences  which  have  been  attributed  to  round  off 
effects.  These  calculations  are  shown  in  Figs.  1  and  2  which  represent 
contour  plots  of  pressure,  density,  energy,  and  speed.  The  time  it  takes 
to  run  this  particular  sample  case  is  30  seconds  on  the  CRAY  and  1.25 
hours  on  the  PC. 

C.  APPLICATION  TO  LINUS 
1.  Initial  Parameters 

Hi  storical  ly [4, 5,6] ,  the  focal  length  of  the  lens  used  in  the 
experiment  has  been  reported  to  be  5.6  to  5.8  cm  and  the  laser  power 
absorbed  is  -10  %.  Recent  measurements  in  a  new  conf iguration[7]  show 
that  the  focal  length  is  4.1  cm  and  the  laser  power  absorbed  is  -40  %. 
Since  the  calculation  shown  here  precedes  the  new  determination,  the 
numbers  used  are  all  based  on  the  earlier  measurement. 

The  dimension  of  the  focal  region  where  the  laser  energy  is  deposited 
is  calculated  from  the  focusing  lens  formula  (see  Section  V.B).  With  a 
5.8  cm  focal  length  lens,  it  results  in  a  focal  radius  of  15  urn  and  length 
of  focal  region  of  110  fxm.  The  measured  initial  laser  energy  is  1  J  but 
^he  use  of  100  mJ  and  greater  energy  input  to  the  code  resulted  in  an 
error  during  runtime  which  still  waits  to  be  resolved.  However,  the  code 
runs  successfully  when  the  energy  input  is  reduced  to  10  mJ.  Therefore, 
the  amount  of  energy  deposited  used  in  the  calculation  is  taken  to  be  10 
mJ  and  is  assumed  to  be  uniform  and  instantaneous.  The  energy  deposited 
is  entered  as  a  multiplier  of  the  initial  specific  internal  energy.  All 
the  units  in  the  code  are  in  the  cgs  system. 
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Fig.  1A.  SALE  Calculation  for  the  Sample  Case 
Running  on  the  CRAY. 
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Fig.  2C,  SALE  Calculation  for  the  Sample  Case 
Running  on  the  PC. 
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Using  the  ideal  gas  law  with  an  initial  pressure  of  150  torr  of 
oxygen  gas  at  room  temperature,  the  initial  mass  density  ' s  2 .4xl0‘‘^g/cm^ . 
The  ideal  gas  equation  of  state  is  assumed  throughout  the  calculation. 

The  computing  mesh  in  cylindrical  coordinates  is  shown  in  Fig.  3. 
This  is  a  20x40  grid  with  each  cell  defined  by  a  10  nm  square.  The  laser 
propagation  axis  is  the  left  boundary  which  is  also  the  symmetry  axis  in 
the  cylindrical  coordinates.  Laser  energy  of  10  mJ  is  deposited  into  a 
focal  region  of  10  ^m  by  120  /im.  An  implicit  Eulerian  calculation  is 
performed  up  to  5  nsec. 

2.  Results 

Figs.  4,  5,  and  6  show  the  results  of  the  calculation  at  1,  3,  and 
5  nsec  respectively.  One  can  clearly  see  the  plasma  expansion  and  the 
shock  front  propagation  as  time  proceeds.  The  initial  region  for  energy 
deposited  is  cigar-shaped  and  it  quickly  expands  spherically  to  0.4  mm  in 
less  than  5  nsec,  thus  the  full  size  of  the  spark  (2  mm)  may  be  reached  in 
-20  nsec.  This  implies  that  it  may  not  be  feasible  to  use  the  optical 
multi-channel  analyzer  (OMA),  which  has  a  30  nsec  time  resolution,  to 
study  the  expansion  behavior  of  the  plasma. 

The  density  plots  show  a  drop  of  about  an  order  of  magnitude  in  the 
core  region.  This  means  that  the  effect  of  the  laser  energy  deposition  is 
to  rapidly  sweep  the  gas  outward,  as  one  intuitively  expects.  In 
addition,  the  figures  show  that  the  shock  front  velocity  is  on  the  order 
of  10®  cm/sec.  This  may  be  compared  to  values  obtained  by  other 
researchers.  Wilke[8]  reported  a  value  of  5.7x10®  cm/sec  in  50  torr  air. 
Askar'yan[9]  gave  a  value  of  10®  cm/sec.  Laser/foil  target  plasma 
experiments  in  NRL[10]  measured  a  value  of  3x10^  cm/sec  in  15  mTorr  H2 
with  10  J  of  laser  energy,  although  this  experiment  was  with  a  solid 
target  implying  different  energy  accomodation.  The  value  presently 
calculated  is  thus  very  close  to  previously  reported  numbers.  The 
velocity  of  sound  under  the  initial  conditions  is  about  4x10^  cm/sec. 

Blast  wave  theory[ll]  gives  a  formula  which  relates  the  shock  front 
velocity  to  the  pressure  behind  shock,  which  is 
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Fig.  3  The  Computing  Mesh  being  Input  into  the  SALE  Code. 
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4C.  SALE  Calculation  for  the  LINUS  Case  at  1  nsec. 
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where  =  pressure  behind  shock 
7  =  ratio  of  specific  heats 
Pq  =  initial  density 
r  =  shock  front  radius 
t  =  time. 

Calculation  of  the  shocK  front  velocity  at  various  times  using  the  above 
equation  is  shown  in  Table  1. 


TABLE  1 

6  6 

T{nsec)  V(shock) (SALE)  10  cm/sec  v(shock) (theory)  10  cm/sec 
1  3.5  3.6 

3  1.8  1.9 

5  1.2  1.4 


We  can  clearly  see  the  excellent  agreement  between  the  shock  velocity  from 
SALE  and  from  theory  --  the  difference  is  less  than  20%.  The  reason  for 
the  close  agreement  is  that  modified  blast  wave  theory  does  not  have  any 
loss  mechanism  and  it  assumes  a  constant  ratio  of  specific  heats  -- 
assumptions  which  are  also  used  in  the  present  SALE  calculation,  but  which 
will  be  modified  for  a  more  realistic  treatment. 

The  calculations  at  5  nsec  show  distinctly  the  reflection  due  to  the 
presence  of  the  boundary.  We  plan  to  couple  the  code  to  a  regrid  routine 
and  perform  a  regrid  with  coarser  grid  size  before  the  shock  wave  reaches 
the  boundary  in  order  to  avoid  any  boundary  effects.  Another  method  would 
^0  start  with  coarse  grid  size,  but  the  limit  on  the  coarsest  grid  size 
is  the  focal  radius  of  the  laser  beam,  i.e.  15  fj,m. 

From  the  pressure  and  density  output,  one  calculates  the  kinetic 
temperature  of  heavy  particles  behind  the  shock  front  to  be  10°  K.  One 
can  also  calmlatp  this  temperature  using  the  velocity  plots  since 
velocity  is  proportional  to  the  square  root  of  the  temperature.  Thi^ 
value  is  an  upper  limit  to  the  temperature  since  no  radiation  lov^. 
processes  have  yet  been  included. 
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3.  Spark  Dimension 

In  order  to  obtain  an  estimate  of  the  time  scale  of  emitted 
radiation,  Sowle  and  Sappenf ield[12]  have  suggested  a  crude  calculation. 
This  results  in  an  estimate  of  the  spark  dimension  which  is  very  close  to 
the  actual  observed  spark  size.  The  assumptions  made  in  this  calculation 
are:  (1)  all  excitations  result  in  radiative  transitions,  (2)  no 
distinction  is  made  among  the  different  ionic  states  of  atomic  oxygen,  and 
(3)  the  emitted  photon  energy  is  10  eV.  This  means  that  the  oxygen  atoms 
are  emitting  radiation  as  hard  as  possible,  i.e.  one  photon  is  emitted  per 
col  1 i b i 0  n . 

Since  the  largest  possible  electron  excitation  rate  coefficient 
expected  is  -10'®  cm®/sec  and  the  assumed  photon  energy  is  10  eV,  the 
emission  coefficient  is  -10'^  eV  cm®/sec.  Multiplying  by  the  oxygen 
density  of  5x10^®  cm'®  (P=150  torr  at  room  temperature)  gives  an  energy 
emission  rate  of  1.5x10^®  eV/sec.  The  SALE  calculation  is  performed  up  to 
5  nsec,  by  which  time  the  energy  expended  is  SxlO'^eV. 

The  energy  being  input  to  the  gas  is  400  mJ  (using  the  new 
measurement).  Energy  deposition  calculation  (see  Section  V.C)  shows  that 
this  results  in  an  energy  input  of  1.3x10^  eV/atom.  The  SALE  calculation 
for  150  torr  pressure  shows  that  the  shock  front  radius,  r,  at  5  nsec  is 
200  We  know  that  the  energy  loss,  which  is  a  measure  of  volumetric 
radiance,  is  directly  proportional  to  r® .  Therefore,  if  we  were  to 
increase  the  radius  by  a  factor  of  10,  i.e.  to  r=2  mm,  the  energy  loss 
would  increase  by  a  factor  of  1000.  At  this  radius,  the  energy  loss  is 
approximately  equal  to  the  energy  input.  Hence,  the  calculation  shows 
^hat  substantial  energy  dumping  would  not  occur  until  the  shock  radius  is 
2  mm.  This  value  is  reasonably  dose  to  the  observed  spark  size.  Thus  it 
appears  from  this  calculatiou  that  the  assumption  of  no  energy  loss  a’ 
oar  1 y  times  is  valid. 

D.  RECOMMENDATIONS 

Important  mfor.mation  was  gained  in  this  study  of  tho  discharqo 
region.  However,  this  initial  study  has  not  accounted  for  many  processes 
which  are  known  to  occur  in  the  LINUS  experiment,  of  which  the  main  ones 
ore  the  optical  loss  processes  and  the  presence  of  plasma  chemistry. 
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Radiation  loss  is  a  major  cooling  channel  in  an  actual  airburst  and 
accounts  for  about  80%  of  the  energy  loss[13].  In  the  laser  produced 
plasma,  uv  and  x-ray  emissions  are  major  cooling  mechanisms  which  are 
certain  to  reduce  the  temperature  behind  the  shock  front.  The  SALE  code 
in  its  present  form  does  not  allow  for  any  energy  loss.  For  it  to  be  more 
applicable  to  the  analysis  of  the  LINUS  experiment,  it  is  necessary  to 
account  for  energy  loss  due  to  dissociation,  exci tat i on/ i on i zat i on ,  and 
chemical  reactions.  In  addition,  we  plan  to  include  uv  and  x-ray 
deposition  and  pro  dosing  of  the  gas  through  which  the  shock  propagates. 

Laser  induced  plasma  contains  more  than  one  species.  For  the  LINUS 
experiment  we  have  a  pool  of  electrons,  positive  atomic  ions(the  highest 
charged  species  identified  is  and  oxygen  i  it 'al  .  A  tremendous 

amount  of  chemistry  goes  on  among  these  species.  Ail  these  chemical 
species  should  be  incorporated  into  the  code.  In  addition,  the  ideal  gas 
equation  of  state  does  not  account  for  i  ntermol ecul ar  force^  and  is 
-triclly  only  valid  at  low  pressure  and  temperature.  Since  the  I.  INUS 
plasma  is  high  in  both  pressure  and  temperature,  a  more  realistic  equation 
-if  state  should  oe  incorporated. 

Another  feature  which  has  been  disregarded  in  the  calculation  is  the 
laser  pulse  shape  and  duration.  The  present  calculation  assumes  an 

instantaneous  deposition.  However,  the  result  shows  that  the  shock  is 
expanding  so  rapidly  toward  the  laser  that  the  high  density  shock  front 
might  interact  further  with  photons  present  at  later  time.  Previously[4] 
the  pulse  shape  of  the  laser  has  been  measured  to  be  spikey  in  normal 

operation  with  no  etalon  in  the  oscillator  cavity.  The  pulse  width  given 
p specification  is  10  nsec.  Because  the  pulse  shape  is  spikey.  the 

effective  pulse  width  is  most  likely  less  than  10  nsec.  A  new 

:et ecm I nat i on [ 7 ]  gives  a  pulse  width  of  20  nsec  and  a  smooth  profile  with 
d  rm*'  time  of  -2  nsec  and  decay  time  of  -5  nsec.  The  new  measurement  is 
■n  tictrast  tn  the  old  one  and  may  be  due  to  an  unaccounted  time  constant 
•'1  * '  ■'  sy  .tr--n.  It  is  suggested  tha^  the  measurement  be  repeated  and  the 
'ir-'c';  laser  pulse  shape  be  incorporated  into  the  SALE  code. 

A  rngrid  routine  will  be  coupled  to  the  code  to  calculate  expansion 
■.t-  later  times.  Ihe  purpose  of  the  regrid  routine  is  to  minimize  the 
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number  of  cells  in  the  calculation  and  hence  shorten  the  runtime  of  the 
problem.  Initially,  one  can  use  a  small  number  of  cells  with  fine  grid. 
As  the  expansion  develops  and  approaches  the  boundary,  a  regrid  with 
coarser  grid  will  allow  the  expansion  to  continue  with  no  interference 
from  the  boundary. 

To  study  whether  the  structures  observed  on  the  contour  plots  are 
authentic  (see  Figs.  4-6),  one  must  perform  sensitivity  tests  to  see  if 
they  change  as  a  function  of  grid  size  or  number  of  cells  or  if  some  other 
parameters  or  effects  are  responsible  for  them. 

The  shape  of  the  spark  has  been  measured  at  200  torr  to  show  evidence 
for  a  donut  around  the  waist  of  the  discharge[5] .  Since  the  focal  region 
of  the  laser  is  bow  tie  shaped,  the  donut  could  arise  due  to  interaction 
of  shock  waves  moving  away  from  the  focal  region  (see  Fig.  7).  It  is 
recommended  that  the  rectangular  focal  region  input  to  the  SALE  code  be 
replaced  by  a  "bow  tie"  shaped  region,  which  is  closer  to  reality,  in 
order  to  test  for  effects  of  shock  wave  interaction.  Should  this 
interaction  occur  in  the  calculation,  it  suggests  a  major  set  of 
relatively  inexpensive  experiments  for  LINUS  in  investigations  of 
interacting  shock  fronts. 


EXPANDED  VIEW  OF  FOCAL  REGION 


Fig.  7.  Bow  Tie  Geometry  for  Laser  Focal  Region  Illustrating 
Shock  Wave  Interaction. 
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III.  STARK  BROADENING  IN  HIGH  ANGULAR  MOMENTUM  STATES  OF  ATOMIC  OXYGEN 
A.  INTRODUCTION 

This  section  describes  research  leading  to  an  understanding  of  line 
broadening  of  high  angular  momentum  states  of  atomic  oxygen  emitted 
radiation  under  plasma  conditions.  The  experiments,  based  on  the  LINUS 
experiment,  yield  time  resolved  spectra  in  the  infrared  in  the  4  ixm  to 
10  iim  spectral  region.  It  is  hoped  that  results  will  lead  to  measurements 
of  electron  densities  of  the  plasma  electrons,  and  an  understanding  of 
properties  of  high  angular  momentum  states  of  atomic  oxygen  and  its  ions. 
For  this  study  we  have  computed  dipole  matrix  elements  for  an  expanded 
basis  set  of  up  to  18  oxygen  atom  energy  levels  for  the  quintet  states 
(6h;  6g;  6f;  5g;  5f;  4f;  7i;  7h;  7g;  7f;  5d;  4d;  8f:  8g;  8h ;  8i  ;  8 j ;  6d), 
and  have  worked  on  an  improved  lineshape  program. 

The  lineshape  function  is  given  by 

S(.)  •  4-  I 

where  in  second  order 

,  ,  47fn_  /-O')  o  ,  f  f(v)  r  dp 

‘I'aa  =  -  — [-J2^<a|r^|a"><a"|r^|a'>  dv  -l[A(z,z'  +  iB(z,z' )  ] , 

and  where  [A(z,z'  +  iB(z,z')]  represents  electron  atom  collision 

integrals,  f(v)  the  electron  velocity  distribution,  W(e)  the  electron 
microfield  distribution,  and  other  terms  have  their  usual  meanings.  The 
term  represents  a  linewidth  term. 

We  are  working  to  apply  this  expression  to  the  near  degenerate  cases 
of  high  principal  quantum  number  and  high  angular  momentum  states.  At  the 
present  time,  electron  collision  processes  are  represented  by  a  single 


''elaxation  time,  T^  and  the  observed  linewidth  is  caused  by  this 
relaxation  time,  the  spontaneous  emission  time,  and  instrument  resolution. 
Because  of  the  near  degeneracy  and  the  great  number  of  J  substates,  and 
therefore  the  sensitivity  of  the  system  to  electric  field  perturbations. 
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our  spectral  model  consists,  at  this  stage,  of  the  large  number  of 
transitions  possible  from  the  perturbations  of  the  electric  microfield  on 
these  high  n,  high  1  levels. 

The  interaction  Hamiltonian  for  the  atom  in  the  presence  of  a  static 
field  is  given  by 

H  =  Hq  +  ezE 

where  Hq  is  the  Hamiltonian  for  the  free  atom.  The  eigenvalues  for  these 
states  are  the  known  energy  levels  of  atomic  oxygen  while  the 
eigenfunctions  are  unknown.  We  have  assumed  hydrogenic  type 
eigenfunctions  as  a  first  approximation.  The  states  under  consideration 
yield  a  number  of  transitions  in  the  infrared  from  2  fim  Lo  16  jum.  Fig.  8 
shows  these  transitions  for  zero  plasma  field  and  18  basis  states. 
Transition  probabilities  are  obtained  using  appropriate  hydrogenic  type 
wave  functions.  In  the  presence  of  the  electric  microfield  these  basis 
states  are  connected  by  dipole  matrix  elements  and  the  states  mixed,  or 
scrambled,  becoming  linear  combinations  of  each  other,  |'l'^>=2Q|(/iQ><n|  i> 
where  <nvi>  are  mixing  coefficients  relating  the  eigenfunctions  of  the 
Stark  problem  to  the  basis  states.  These  eigenfunctions  are  determined  by 
diagonalizing  the  Stark  matrix.  Because  of  state  mixing,  previously 
forbidden  transitions  are  now  possible.  The  strengths  of  these  forbidden 
transitions  will  depend  on  the  Stark  field  through  the  mixing 
coefficients.  For  example,  the  transition  "6h"  -  "5f"  is  proportional  to 

where  the  6h,  5f  refer  to  parent  states  that  are  now  "mixed"  by  the  plasma 
electric  field  with  other  states  and  where  the  matrix  elements  for  the 
basis  states  have  been  calculated  as  a  part  of  the  Stark  effect 
interaction.  The  mixing  coefficients,  as  mentioned,  are  field  dependent. 
The  work  done  has  been  concerned  with  the  calculation  of  these  dipole 
matrix  elements  between  all  appropriate  JMj  states,  the  extention  of  the 
calculation  to  larger  basis  sets,  the  inclusion  of  J  substates,  and 
finally  the  inclusion  of  Mj  substates. 

32 


microns 


6  States 


Zero  Plasma  Field  Spectra  for  18  Basis  States  -  6-16  microns 


w-  M-  V-  -  -  .  -y-if-k—  I.  r-  ^.■.  ■■■■X-  V  «  -  a  »1l  »  J  ■  v»1 


B.  STARK  CALCULATION 

We  have  divided  the  Stark  problem  into  two  parts,  the  electron 
microfield  that  the  atoms  interact  with,  and  the  resulting  atomic  states 
responsible  for  transitions. 

1.  Atomic  State  Matrix  Elements 

We  must  calculate  all  matrix  elements  <0n'J'M'L'sl'^M^‘^nJMLS^ 
using  appropriate  wavefunctions.  These  matrix  elements  will  be  used  to 
construct  the  matrix  for  the  Stark  Hamiltonian, 

^o^^'^n'J'M'L'sl  V'^nJMLS^ 

and  they  will  also  be  used  to  compute  radiative  transition  probabilities. 
Diagonal ization  of  this  matrix  results  in  eigenvalues  and  eigenfunctions 
for  a  given  electric  field.  Thus,  for  each  field  we  obtain  eigenenergies 
and  transition  probabilities  and  thus  a  Stark  spectrum.  The  convolution 
of  all  these  thousands  of  lines  using  the  fundamental  widths  caused  by  Tj 
processes  results  in  a  Stark  broadening  line  profile. 

2.  Tne  Electric  Microfield 

In  order  to  find  the  electric  field,  e,  we  need  to  consider  the 
problem  of  a  fluctuating  charge  density  (the  plasma  in  our  case).  There 
are  two  traditional  approximations  where  either  the  ions  are  stationary 
relative  to  the  electron  motion  (high  frequency  approximation)  or  the  ion 
motion  dictates  the  time  variation  of  the  system  (low  frequency 
approximation).  The  essential  physical  picture  is  that  of  a  sphere  whose 
radius  is  given  by  the  Debye  length,  =  (kTg/47rnge^)^'^^,  the  "thickness" 
of  the  ionic  atmosphere  surrounding  the  point  of  interest  (centered  at  the 
emitter)  through  which  electrons  travel.  The  electrons  of  density  ng  have 
a  mean  radius  r^  =  [3/47rng]  .  An  electron  penetrating  the  Debye  sphere 

at  a  distance  Cg  from  the  emitter  provides  an  electric  field  of  strength 
e/rg^.  The  volume  containing  such  charges  (and  hence  the  number 
distribution)  depends  on  the  radius  rg.  Electrons  outside  the 

characteristic  distance  are  exactly  balanced  by  the  background  of  positive 
charges.  The  actual  condition  is  quite  complicated  so  we  shall  rely  on 
the  computations  of  Baranger,  Hooper  and  others  for  values  of  the  field 
distribution.  The  parameter  of  the  distribution  is  the  ratio  of  the  mean 
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electron  distance  to  the  Debye  length,  a  =  function  of  the 
electron  number  density,  n^,  and  the  electron  temperature,  Tg)  and  the 

O 

reduced  electric  field  variable  e/e„,  where  ^  eA',  .  We  have  chosen  to 
fit  the  points  calcu''a'ed  by  Hooper  to  a  convenient  function  for  use  in 
the  Stark  calculation.  The  Stark  calculation  begins  by  choosing  the 
electron  number  density  n^,  the  electron  temperature  Tg,  an  electronic,  or 
spectroscopic,  temperature  Ig-j  with  subsequent  calculation  of  the 
microfield  distribution,  computation  of  the  Stark  matrix  followed  by 
siatnix  d i agonal i zat i on .  For  each  point  on  the  microfield  distribution 
there  corresponds  a  probability  of  occurrence  and  a  reduced  electric  field 
strength.  The  field  at  the  emitter  is  the  reduced  field  multiplied  by  the 

p 

mean  field  fg  =  field  becomes  the  I'-i'id  part  of  the  Stark 
matrix.  Subsequent  di agonal ization  gives  eigenvalues,  transitions  and 
transition  probabilities 

R  =  W(Q;,€/fQ)  i<’F^ 

for  eacn  point  of  the  microfield  distribution.  These  data  are  introduced 
into  a  line  shape  program  that  assumes  a  Lorentzian  shape  whose  area  is 
proportional  to  the  transition  probability  weighted  by  an  appropriate 
Boltzmann  factor  and  whose  width  is  given  by  relaxation  processes  and 
instruiTient  resolution  (LT  =  l/T^  r  NgOV  +  ...). 

The  results  for  18  states  are  shown  in  Fig.  9  for  several  electron 
densities.  These  are  preliminary  results  since  a  number  of  approximations 
and  omissions  exist.  First,  we  are  considering  only  S  =  2,  the  quintet 
states;  second,  we  hdve  not  taken  degeneracies  into  account,  nor  in  these 
spectra,  existence  of  i  sublevels.  Thus,  the  intensities  are  relative. 
Fur’-.riermore ,  we  have  parameterized  electron  collisions  treating  the 
problem  a.  a  static  one  consisting  of  average  fields  at  the  emitters.  The 
essential  conclusion  is  that  the  line  positions  are  very  sensitive  to 
electron  microfield  strength,  as  expected  from  the  near  degeneracy  of  many 
of  tne  energy  levels.  From  a  comparison  of  these  results  with  the  LINUS 
aata,  we  conclude  that  the  electron  density  is  in  the  neighborhood  of  10^'^ 
el ectrons /cm^  at  a  delay  time  of  6  /iS  after  initiation  of  the  plasma. 
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IV.  KINETIC  EFFECTS 

Preliminary  modelling  of  the  kinetic  effects  has  been  described  in 
detail  [4],  The  initial  concentrations  assumed,  based  on  energetics,  are 
80%  0'*',  15%  and  5%  0^'*’.  The  reaction  set,  which  includes 

recombination  and  reionization,  is  shown  in  Table  Z.  The  12  kinetic 
equations  are  solved  by  PCCHEM,  a  numerical  integration  routine  written  by 
Gene  Adams  at  Utah  State  University  which  has  been  modified  by  Armstrong 
for  this  calculation. 

Previously[5] ,  the  time  history  of  the  species  concentration  was 
calculated  by  assuming  an  exponential  temperature  profile.  We  have  now 
performed  a  new  calculation  using  a  new  temperature  profile  obtained  semi- 
empirically  by  assuming  that  the  electrons  are  in  collisional  equilibrium 
with  the  dominant  charged  state.  This  means  that  the  electrons  are 
coupled  with  the  highest  charged  state  observed  at  early  time  (i.e. 
and  with  the  lowest  charged  state  observed  (i.e.  O'*”)  at  late  time.  The 
temperature  profile  is  then  determined  from  the  latest  available  data  on 
spectroscopic  temperature[14]  shown  in  Table  3.  From  examining  this 
table,  one  observes  that  at  early  time  the  O^'*’  temperature  is  on  the  order 
of  10^  K  and  decreases  sharply.  At  100  nsec  the  O'*'  temperature  is  about 
3x10^  K  and  from  then  on  varies  slowly.  The  temperature  for  0  neutral  at 
late  time  is  determined  to  be  -8000  K[15].  Thus,  we  assumed  a  temperature 
profile  which  has  an  exponential  function  from  0  to  100  nsec  and  behaves 
linearly  from  100  nsec  to  1  /xsec.  The  analytic  form  of  the  new 
temperature  profile  is 


T  =  TQ[exp(-t/2.7xl0^)-2.2xl0^t+0.3] 


and  is  shown  in  Fig.  10. 

The  results  of  the  new  calculation  at  10  and  100  torr  are  shown  in 
Figs,  11  and  12  respectively.  The  general  trend  of  the  time  profile  of 
the  species  concentration  is  the  same  as  before.  One  can  see  that  0 
neutral  does  not  dominate  until  t>0.2  /isec  at  which  time  0^’*'  and  0^"^  are 
no  longer  important.  Again  we  see  the  dip  in  0'*'  at  t>l  nsec  accompanied 
by  an  increase  in  0^'*'  which  is  attributed  to  reionization  of  O'*'. 
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TABLE  2 


REACTION  SET  FOR  LINUS  RECOMBINATION 


RXN 

k(T)»A(T/300)‘’ 

e-c/T 

1. 

O'^'^  +  e  +  e  --> 

0^+  +  6 

1.9E-18 

-4.5 

0 

2. 

0^"*"  +  e  --> 

02+ 

*^d 

6.2E-11 

-1.5 

5.6E3 

3. 

0^"^  +  a 

0^+ 

3.2E-11 

-0.7 

0 

4. 

0^"*"  +  e  +  e  --> 

O'*’  +  e 

5.6E-19 

-4.5 

0 

5. 

0^"^  +  e 

O'*’ 

•^d 

2.8E-11 

-1.5 

5.6E3 

6. 

0^^  +  e  --> 

O'*" 

1.4E-11 

-0.7 

0 

7. 

0"^  +  e  +  e  --> 

0  +  e 

7.0E-20 

-4.5 

0 

8. 

O'*"  +  e  --> 

0 

^d 

6.9E-12 

-1.5 

5.6E3 

9. 

O'*"  +  e  --> 

0 

3.5E-12 

-0.7 

0 

10. 

0  +  e  --> 

0'*’  +  e  +  e 

*‘‘1 

1.7E-12 

1.5 

1.6E5 

11. 

0"^  +  e  --> 

0^'*’+  e  +  e 

*^i 

5.4E-13 

1.5 

4.0E5 

12 

0^"^+  e  --> 

0^'*’+  e  +  e 

4.1E-13 

1.5 

6.4E5 

T  -  T, 

[exp(-t/2. 

17E-8) 

-  2.2E5t  + 

0.3] 

ic^athree-body,  k^adielectronic,  k^=radiative,  lc^»re-ioni2ati( 


Because  the  new  profile  has  lower  temperatures  throughout  most  of  the 
time  domain,  two  differences  in  the  new  calculation  are  observed.  First, 
we  see  a  general  trend  of  compression  for  all  species  at  early  time. 
Second,  the  electron  and  O'*'  densities  are  higher  by  a  factor  of  5  at  1 
lisec  in  the  new  result.  The  two  differences  can  be  explained  by  the  fact 
that  the  dominant  three-body  recombination  rate  is  an  inverse  function  of 
temperature.  The  recombination  rate  using  the  new  profile  is  higher  and 
hence  accelerates  the  early  time  build  up  of  lower  charged  species. 
Similarly,  the  temperature  of  the  new  calculation  at  1  Isec  is  higher 
(8000  K  vs.  4000  K) .  Hence,  the  slower  recombination  rate  results  in  a 
higher  electron  and  0'*'  densities  at  late  time. 

The  assumptions  used  in  this  calculation  may  be  incorrect  but  they 
give  us  a  good  place  to  begin  the  model  and  will  help  us  in  pursuing 
further  more  detailed  modelling  which  will  include  time  dependent  terms 
for  all  the  excited  states  in  each  species. 
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V.  MISCELLANEOUS  ISSUES 

A.  Spherical  Aberrations 

The  purpose  of  this  calculation  is  to  determine  whether  the  lens  used 
in  the  LINUS  experiment  is  a  source  of  spherical  aberration. 

The  number  of  wave’jngths  of  spherical  aberration  is  given  by[16] 

'l'(p)  -(ap)^^  n  (n+1)^ 

m  =  - = - -  [ - - - +  - ] 

A  8Af-^  4(n-l)^  4(n+2)  n(n+2) 

where  m  =  aberration  number(wavelengths) 

p  =  r/a  =  radial  distance(cm)  from  optical  axis  at  front  lens 
surface  divided  by  lens  radius(cm) 

A  =  wavel ength (cm) 

n  =  index  of  refraction  of  lens  material 
f  =  lens  focal  length{cm). 

For  n  =  1.55,  this  reduces  to 
m  =  §(r)/A  =  -3.056rV8Af^. 

For  r^3^  =  beam  radius  =  0.3  cm, 

m  =  -0.166  wavelengths  (minus  sign  means  towards  the  beam  source). 

This  result  is  consistent  with  figure  9.2,  p.  722  from  Morgan's 
article[16]  which  shows  that  for  a  6  mm  beam  at  1.06  wavelength,  '1><0.2 
at  a  5.6  cm  focal  length.  For  a  4.1  cm  focal  length  lens,  m  =  -0.42 
wavelengths.  Thus  the  calculation  shows  that  spherical  aberration  caused 
by  the  present  lens  is  negligible. 

B.  Focal  Volume  Description 

The  radius  of  the  focal  volume  is  given  by  the  thin  lens 
approximation  from  geometrical  optics[17].  The  thin  lens  approximation  is 
valid  if  the  lens  thickness  is  small  with  respect  to  the  focal  length. 
For  the  LINUS  lens,  this  ratio  is  <0.1  so  the  thin  lens  approximation  is 
probably  valid.  In  this  case  one  obtains  for  the  beam  diameter 


tan^/2  =  (d/2)/f 


f.' 

f.' 

K' 


i 

y. 

y. 

N‘. 

y.' 
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where  f  =  lens  focal  length 

6'2  =  half-angle  laser  beam  divergence(radians) 
d '2  =  hal f-di ameter  focal  waist 
which,  for  small  angles  reduces  to 
d/2  =  r  =  f{8/2) 

For  LINUS,  9/2  ~  0.25  milliradians,  and  f  =  5.8  cm,  thus  r  =  15  /xm. 

For  the  focal  length,  1,  the  general  definition  is  given  as  that 

distance  where  the  intensity  drops  to  half  of  that  at  the  focal  plane. 

Again,  from  the  thin  lens  approximation  and  geometrical  optics 
arguments [  16] , 

1  =  2(72  -  Uf^^/D 

where  the  variables  are  defined  above  and 

0  =  unfocussed  beam  diameter  at  the  lens  surface. 

For  our  case,  D  =  6.5  mm,  thus  1  =  110  /xm. 

For  the  new  focal  length  measurement,  f  =  4.1  cm,  r  =  10  /xm  and 
1  =  54  ^m. 

These  numbers  are  valid  to  the  extent  that  the  thin  lens 

approximation  is  appropriate  and  to  the  extent  that  self-focussing  is 
unimportant.  However,  any  uncertainties  in  these  numbers  are  likely  to  be 
much  less  than  uncertainties  in  other  areas  that  are  more  critical  to  our 
analysis . 

C.  Energy  Deposited  per  Molecule 

Given  an  initial  laser  power,  it  is  instructive  to  calculate  the 
maximum  amount  of  energy  which  could  be  deposited  onto  each  O2  molecule. 
The  number  of  molecules  in  the  laser  focal  volume,  n,  is  given  by 

n  =  NV, 

where  N  =  number  density  of  molecules,  and 
V  =  focal  volume  of  laser. 

For  a  pressure  of  150  torr  at  room  temperature,  N=4.8xl0  cm  The 

focal  volume,  V,  is 
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V  =  Ttr^l  . 

whet'e  r  =  focal  t'ddiiis,  and 

1  =  length  of  the  focal  region. 

The  two  parameters,  r  and  1.  can  be  calculated  using  the  focusing  lens 
formula  (see  Section  V.B).  For  the  lens  used  in  the  experiment  (focal 
length  of  4.1  cm),  r=10  iim  and  1=54  /im.  This  results  in  a  focal  volume  of 
2x10  cm  .  Flence,  the  number  of  O2  molecules  in  the  focal  volume  is 
-10^^.  The  initial  laser  energy  is  ~i  J  per  pulse  and  typically  40%  of 
this  energy  is  absorbed[7].  Therefore,  the  maximum  energy  deposited  is 
-10^  eV  per  O2  molecule. 

The  dissociation  energy  is  5.2  eV.  The  first,  second,  and  third 
ionization  potential  of  0  atom  are  respectively  13.6,  35. I*,  and  55.0  eV. 
Thus  the  generation  of  two  0^'*'  ions  from  one  O2  molecule  requires  energy 
of  213  eV.  This  quantity  is  a  minute  fraction  of  the  initial  laser 
energy.  This  implies  that  most  of  the  energy  is  going  into  PV(pressure 
volume)  work. 

D.  Momentum  Transfer 

Another  issiie  which  needs  to  be  addressed  is  the  shift  of  the  spark 
location  as  a  result  of  momentum  transferred  from  the  high  power  laser 
beam  to  the  gas.  The  momentum  of  a  photon,  M,  is  given  by 


M  =  2E/V, 

where  £  =  hu  is  tlie  photon  energy  and  v  =  c  is  the  velocity  of  light.  For 
photon  of  wavelength  1,06  /urn,  M=1.3xl0"^^  g-cm/sec. 

The  distance  shifted  due  to  photon  momentum  imparted  to  the 
mo''ecules,  ax,  is 


AX  =  MAt/m, 

where  At=effective  pulse  width  is  assumed  to  be  2  nsec  and  m-mass  of  0; 
moiecules  in  the  focal  volume(P=150  torr  and  T-298  K)  is  5x10  g. 


Therefore,  ax  is  5.2x10'^^  cm  per  photon.  Since  the  amount  of  energy 
contained  in  a  1.06  iim  phoion  is  2x10'^^  erg,  the  initial  laser  energy  of 
400  mj  is  equivalent  to  2x10^®  photons  per  pulse.  The  final  distance 
shifted  due  to  momentum  transfer  should  be  1  mm;  however,  no  such  shift 
has  been  detected  by  eye.  In  fact,  the  shift  may  be  offset  by  laser-shock 
front  interactions  which  cause  the  spark  to  elongate  toward  the  laser. 
Thus  this  issue  is  yet  not  fully  resolved. 

E.  Continuum 

If  data  on  early  time  continuum  can  be  described  by  bremsstrahl ung  in 
the  field  of  an  ion,  then  values  for  electron  density  and  temperature  can 
be  obtained.  The  spectral  emission  coefficient  for  free-free  transitions 
in  the  field  of  an  ion,  is  given  by[ll] 


Ng  exp(-h;//kT)  ^  Z., 


2  + 


where  m  =  electron  mass 

k  =  Boltzmann's  constant 
T  =  electron  gas  temperature 
e  =  electron  charge 
c  =  velocity  of  1  ight 
Ng  =  electron  density 
h  =  Planck's  constant 
1/  =  frequency  of  transition 
=  charge  of  ith  ion 
N:'''  =  density  of  ith  ion. 

3+ 

Since  the  highest  charged  state  observed  in  the  line  spectra  is  0  ,  we 

assume  that  early  time  continuum  is  also  dominated  by  With  this 

assumption,  the  summation  in  the  above  is  reduced  to  one  term:  ZpS  and 
N^‘*'=(l/3)Ng.  The  equation  is  then  simplified  to 


=  constant  x  Ng 


2  T-l/Z 


exp(  -hi'/kT) 
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Fig.  13  is  a  plot  of  the  emission  coefficient  on  a  relative  scale  as  a 
function  of  wavelength.  The  spectral  region  covered  is  3000  to  5000  A 
with  the  electron  density  being  10^^  cm'^.  We  can  see  from  the  figure 
that  the  functional  form  of  the  emission  coefficient  changes  for  different 
values  of  the  electron  temperature.  Therefore,  by  obtaining  emission 
coefficients  from  experimental  continuum  spectra,  we  can  compare  the 
functional  form  of  the  experimental  coefficients  to  sets  of  curves  of 
calculated  coefficients  for  different  temperatures.  The  calculated  curve 
which  most  closely  matches  the  observed  form  will  give  value  for  the 
electron  temperature. 

For  a  fixed  temperature  the  effect  of  varying  the  electron  density  is 
to  change  the  absolute  scale  of  the  emission  rot-.ff  i  ci  ent  with  the 
functional  form  remaining  unchanged.  Thus  to  obtain  elec'^^ron  density  one 
would  need  an  absolute  calibration  of  the  experimental  emission 
coefficient. 

F.  CVF  vs.  SPEX 

The  purpose  of  this  calculation  is  to  determine  whether  there  is  any 
real  advantage  in  using  the  OCLI  CVF  over  the  SPEX  1870  spectrometer  in 
the  LINUS  experiment.  Two  approaches  are  used:  one  calculates  the  net 
radiance  reaching  the  detector  and  the  other  calculates  the  net 
throughput . 

1.  Net  radiance 
a.  Spectrometer 

The  SPEX  1870  is  a  0.5  m  monochromator  which  consists  of  two  mirrors 
and  a  grating  in  a  Czerny-Turner  configuration.  Its  f/#  is  6.9  and  is 
related  to  the  solid  acceptance  angle,  Q,  by 


Q  =  ^/(4f/fi^)  =  0.017  sterad. 


The  slit  width  used  in  the  LINUS  experiment  is  3  mm  and  this  corresponds 
to  nominal  resolutions  of  0.02  /im  in  the  SWIR  and  0.04  /urn  in  the  MWIR  and 
LWIR. 

If  we  start  with  a  point  source  of  radiance,  Np.  and  using  a  lens 
which  matches  the  f/^  of  the  specti'ometer ,  the  radiance  entering  the 
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monochromator  is  (n/47r)Np  =  (l/16f/#^)Np.  Assuming  that  the  loss  of  the 
grating  is  18-25%  and  the  loss  of  each  mirror  is  5-15%,  the  net  radiance 
reaching  the  detector,  N^p,  is 

Nsp  =  Rf;|^R|vi2RG (0/471  )Np,  where  R  =  reflectance, 
best  case  N^p  =  (0. 95)2(0. 82){0/47r)Np  =  9.7xlO'^Np. 
worst  case  N^p  =  (0 .85)2(0. 75)  (0/47r)Np  =  7.1xlO’'^Np. 

Thus,  the  amount  of  radiance  after  passing  through  the  monochromator  is 
-1/1000  of  the  source  radiance. 

b.  CVF 

The  CVF  which  we  are  considering  consists  of  three  80#  segments 
covering  wavelengths  from  2.5  to  14.5  /xm.  The  front  view[18]  is  shown  in 
Fig.  14.  Its  specifications,  which  are  valid  for  normal  incidence  and 
room  temperature,  are 

\  D  ^cvf 


Segment 

Wavelength(^) 

Hal f-Bandwidth(%) 

Peak  Transmittance(%) 

I 

2. 5-4. 5  (SWIR) 

<1.35 

>25 

II 

4. 4-8.0  (MWIR) 

<1.35 

>30 

III 

7.9-14.5(LWIR) 

<1.8 

>30 

If  we  deviate  from  these  conditions,  typical  wavelength  shifts  as  a 
function  of  incident  angle  and  temperature  are  shown  in  Figs.  15  and  16. 
From  the  half-bandwidth,  we  calculate  the  optimum  resolution  to  be 
0.05(SWIR),  0.08(MWIR),  and  0.20(LWIR)  /xm. 

The  optimum  angular  slit  width  in  degrees  of  a  CVF,  4>Qp^,  assuming  a 
collimated  and  normal  incidence  beam,  is[19] 

'^opt  ^  ^**^*^'^/(^max‘^min) 

where  D  is  the  hal f -bandwidth  and  H  is  the  angular  length  of  linearity  in 
degrees.  Using  the  optimum  slit  width  results  in  maximum  energy 
throughput  without  loss  in  resolution.  The  ^opt'^  shown  in  the 
following  table; 


MOUNTED  2-S-14.5U  NARROW  BANDPASS  INFRARED 
CIRCULAR  VARIABLE  FILTER  SETS  WCV.2.5/14.5-L 


I  ; - i27'aiameter 

I  y  KAtA 


i  segment  ii 

180' 


Thickness:  0.180"  nominal 


Fig.  14.  Front  View  of  CVF. 
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'  V. 


I 


II 


III 


X(/i) 

2.5 

3.5 

4.5 
4.4 
6.2 
8.0 
7.9 

11.2 

14.5 


‘^opt(^eg) 

1.35 

1.89 

2.43 
1.32 
1.86 
2.40 
1.72 

2.44 
3.16 


From  the  average  angular  slit  width  of  2.1°,  the  average  linear  slit 
width  is  calculated  (assuming  radius  of  1.6")  to  be  1.5  mm.  From  Fig.  14 
we  can  see  that  the  useable  slit  height  is  0.3".  Thus  the  optimum  useable 
area  of  the  CVF  with  a  rectangular  slit  is  0.11  cm*^. 

Two  optical  systems  may  be  conceived  to  accommodate  the  CVF  in  the 
LINUS  experiment.  The  first  one  is  to  use  a  single  lens  to  focus  the 
light  onto  the  CVF  (see  Fig.  17).  In  order  to  avoid  wavelength  shift  due 
to  off  angle  incident  beam,  a  f/8  or  f/9  system  is  required.  The  net 
radiance  in  this  case  is 


Ncvf  ■  t„f(l/16f/»2)Np 
=  2-3x10  '^Np. 

Thus  this  CVF  arrangement  has  no  advantage  over  the  spectrometer. 

The  second  arrangement  is  to  use  a  fast  lens,  a  collimator,  and  a 
slit  (see  Fig.  18).  In  this  case  the  amount  of  light  reaching  the 
detector  is 

^cvf  ^  '^cvf^c^  ^  (^sl  i  t^^col  1  imator )^p  • 

If  we  place  the  collimator  at  a  position  such  that  its  effective  diameter 
is  approximately  the  same  as  the  slit  height  (7.62  mm),  then  the  ratio, 

'^sl  it/^col  1  imator’  ™  width.  Assuming  a  f/1 

lens  and  R^.  =  85-95%, 
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Fig.  17.  Optical  System  Using  One  Lens  for  Using  CVF  in  LINUS. 


CVF 

I  n 


Fig.  18.  Optical  System  Using  Lens,  Collimator,  and  Slit  for  Using  CVF 
in  LINUS. 


Thus  if  we  use  the  optimum  slit  width,  the  gain  of  the  CVF  over  the 
spectrometer  is  about  a  factor  of  5  with  a  corresponding  resolution  loss 
of  factors  of  2(SWIR),  2{MWIR),  and  5(LWIR).  If  we  use  the  same  slit 
width  as  the  spectrometer,  i.e.  3  mm,  the  gain  of  the  CVF  is  approximately 
an  order  of  magnitude,  but  with  an  even  greater  loss  of  resolution  than 
above . 

2.  Net  Throughput 

The  throughput,  U,  is  defined  as 
U  =  An,  where 

A  =  entrance  aperture  area 
n  =  acceptance  solid  angle. 

a.  Spectrometer 

The  throughput  of  a  grating  spectrometer,  U^p,  is 

U3P  =  slit  area  *  acceptance  solid  angle 
=  hAx(7r/4f/#2) 

where  h  and  ax  are  the  slit  height  and  width  respectively.  With  h=0.2  cm 
and  ax=0.3  cm,  the  throughput  is  0.001  cm^-sterad.  Taking  into  account 
the  loss  on  the  three  optical  elements,  the  net  throughput  is 

Usp(net)  =  (0. 95)2(0. 82){0. 001) 

=  7.3x10'^  cm2-sterad(best  case). 

U3p(net)  =  S.AxlO"'^  cm2-sterad(worst  case). 

b.  CVF 

The  throughput  of  a  CVF  is  the  :,ame  as  that  for  an  interferometer. 
This  is  given  by 

Ucvf  =  27rA/R 

where  R=>/aX  is  the  resolving  power.  The  net  throughput  is 


Ucvf(net)  =  . 


Using  the  optimum  slit  area  of  0.11  cm^,  the  net  throughput  for  the  three 
segments  is 

Segment  U^y|r(net)  (cm^-sterad) 

I  2.4x10'^ 

II  2.8x10-3 

III  3.7x10-3 

The  throughput  of  the  CVF  is  about  a  factor  of  five  higher  than  that  of 

the  spectrometer.  If  we  use  slit  width  of  3  mm,  the  throughput  gain  of 

the  CVF  is  an  order  of  magnitude.  These  results  agreo  with  the  net 

radiance  results.  The  comments  on  the  resolution  loss  also  apply. 

From  the  above  calculations  we  conclude  that  the  CVF  does  have 
throughput  advantage  over  the  SPEX  1870  but  this  occurs  only  with  a 
simultaneous  reduction  in  resolution.  Thus  if  one  requires  more 
throughput  and  if  the  resolution  is  not  a  critical  factor  in  the 
experiment,  then  the  CVF  might  be  considered  if  a  factor  of  5-10  will  be 
important . 

3.  Recommendations 

The  LINUS  experiment  might  be  improved  upon  by  making  the  following 
two  changes; 

a.  Use  a  focusing  lens  which  matches  the  f/#  of  the  monochromator  will 
help  the  light  collection  efficiency  and  permit  the  slit  width  to  be 
reduced(di scussed  in  more  detail  below). 

b.  The  focal  waist  of  the  focusing  lens,  r,  is  diffraction  limited  and 
is  given  by[20] 

r  =  1.22fX/d. 

Assuming  properly  matched  optics,  r  =  70  ^m  at  X  =  8.5  /xm.  At  present,  a 
3  mm  slit  width  is  used  in  the  LINUS  experiment  which  is  much  larger  than 
the  optimum  slit  width  of  140  /im.  Thus  it  might  be  worthwhile  to  close 
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VI.  EXCEDE  and  ARCTIC 

A.  EXCEDE 

We  have  acquired  the  EXCEDE[21]  and  ARCTIC[22,23,24]  codes  in  order 
to  pursue  their  applicability  to  the  LABCEDE  experiment. 

Traditionally,  energy  deposition  problems  have  been  treated  using  the 
continuous  slowing  down  approxiination(CSDA) ,  by  solving  the  Focker- PI anck 
or  Boltzmann  equations,  or  doing  a  Monte  Carlo  calculation.  The  EXCEDE 
code  (the  special  ARCTIC  code  which  was  applied  to  the  EXCEDE  project) 
employs  a  different  theory  of  energy  deposition  to  calculate  energy 
deposition  rate  from  which  the  radial  distribution  of  electron  density  may 
be  derived.  This  code  was  written  by  Tarr[21]  an'  ■.  important  in  its 
potential  application  to  the  LABCEDE  electron  scattering  problem.  The 
code  was  previously  operated  on  the  VAX  780  and  we  have  been  successful  in 
converting  it  to  run  on  the  PC.  Agreement  between  the  results  of  a  test 
case  which  ran  on  both  computers  is  excellent.  The  theory  is  as  follows. 

If  one  starts  with  an  electron  beam  incident  on  a  gas,  the  code 
calculates  the  dynamics  of  electrons  undergoing  n  number  of  collisions 
with  the  scattering  gas.  The  parameters  of  interest  are  position(x), 
pathlength(l ) ,  energy(E),  and  angles  of  velocity  vector  (^,0).  After  n 
collisions  the  ensemble  of  electrons  has  a  spread  in  the  above  parameters 
which  may  be  represented  by  a  probability  distribution  function,  Z.  If  we 
Know  F,  we  can  calculate  quantities  such  as  flux,  4’,  energy  deposition 
rate,  d^/dt,  and  radial  distribution  of  electron  density,  G(R).  These  are 

4>  =  JoXjFdl 
d^/dt  =  jQXjWFdldEdQ 
G(R)  =  ,|(de/dt)dz 

where  J  -  electron  current  density 
X  =  mean  free  path 

W  =  probability  that  electron  of  energy  E  losing  energy  f. 

The  code  assumes  a  form  for  F  and  uses  moments  to  define  the  assumed 
distribution.  The  moments  are  generated  by  differentiating  the 


:‘iaracteri  St  ic  functions  of  the  random  variables.  To  simplify  the 
Jer  i  vat  i  or'  of  the  equations,  the  assumption  of  many  number  of  scattering 
events  was  made.  It  is  this  restriction  which  hinders  the  applicability 
of  the  code  in  its  present  form  directly  to  the  LABCEDf  electron 
scattering  problem  because  the  LABCEDE  experimental  conditions  give  rise 
to  only  a  few  number  of  collisions.  In  order  to  apply  the  theory  to  the 
LABCEDE  problem,  the  moment  equations  must  be  rederived  without  the 
assumption  of  many  scattering  events  and  the  equations  must  be  recoded. 

'^ig.  19  shows  the  output  of  the  EXCEDE  code  for  the  case  of  4.5  keV 
‘^lectr'ons  and  10  collision  events  for  N2  gas.  We  see  that  the  radial 
electron  density  has  dropped  by  two  orders  of  magnitude  as  the  radius 
varies  from  0  to  65  m.  This  is  in  vast  contrast  to  experimentally 
obtained  radial  electron  density  in  Fig.  20  which  shows  that  the  same  drop 
occurs  over  a  distance  of  0.1  m[25].  This  is  not  surprising  from  the 
discussion  above.  To  apply  the  EXCEDE  code  to  the  electron  scattering 
problem,  we  must  rederive  the  set  of  moment  equations  without  making  the 
assumption  of  many  number  of  scatterings. 

B.  ARCTIC 

The  ARCTIC  code  is  a  general  purpose  routine  designed  for  calculating 
production  rates  and  chemistry  when  a  flux  of  auroral  electrons  bombards 
the  atmosphere  under  ambient  and  disturbed  conditions.  The  section  on 
energy  deposition  is  based  on  the  electron  energy  loss  program  written  by 
Peterson  et  al.[26].  The  code  is  composed  of  a  main  driver  program  with 
49  subroutines  which  control  data  input  and  output,  and  calculate  volume 
production  and  emission  rates.  The  code  was  originally  programmed  to  run 
on  the  CDC  7600.  Subsequently  it  was  modified  to  run  on  the  VAX  780. 
However,  the  core  of  the  program  remains  in  Fortran  66.  This  presents 
most  of  the  problems  when  compiling  the  program  on  the  PC(supports  Fortran 
77).  Linktime  errors  a'-ise  from  calling  routines  which  are  not  supported 
on  the  PC.  Runtime  errors  can  be  attributed  to  either  the  architecture  of 
the  PC  or  idiosyncrasies  of  the  fortran  compiler  (IBM  Professional 
Fortran)  . 
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After  making  all  the  necessary  changes,  the  program  is  tested  by 
comparing  to  the  results  of  a  ELIAS  drizzle  case  from  the  VAX.  Agreement 
between  the  calculations  using  the  two  computers  is  excellent. 
Differences  occur  in  the  last  decimal  place,  i.e.  less  than  1%  error.  On 
the  VAX  this  problen,  typically  takes  ~2  hours  to  complete.  On  the  XT,  the 
production  rates  calculation  takes  -21  hours  and  the  chemistry  code  takes 
1/2  to  1  hour  per  print  time.  On  the  COMPAQ  286,  the  running  time  reduces 
to  17  hours.  The  size  of  the  ex-^cutable  version  of  the  code  is  -600 
Kbytes  which  is  just  oelow  the  PC  memory  capacity  of  640  K. 
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